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Abstract 

Using the representation of numbers, the methods of decreasing the number 
of calculation steps for Pell's equation are developed. The formulae relating the 
natural solution xq of Pell's equation for with the infinite number of minimal 
natural solutions Xi for other Ai are obtained for various cases. The parametric 
representation of Pell's equation solutions are obtained with the help of new 
"method of hidden parameters". 



1 INTRODUCTION 

Pell's equation for natural numbers A, x, y (actually put forward by Fermat in 1657) 
(1) -Ax^ = l 

has been thoroughly studied (see [1-6,9] and references herein). In the subsequent 
text, we shall consider the minimal nontrivial natural solutions {x,y) ^ (0,1) only 
(each solution is unique; the value of x determines the value of y\ To solve (1) for 
the given natural number A which is not a square, one can use several methods: 
the sequential differences (in essence, Euclidian algorithm: measurement of a greater 
number by a smaller one), the method of Wallis-Brouncker, Indian cyclic method, the 
method of continued fraction. We shall not consider the method of Wallis-Brouncker 
and Indian cyclic method, since they are rather unwieldy for using " a pen and a paper 
only" , and they use some exhaustive search: the number of algorithmic steps is great 
enough (the cyclic method can be connected with the method of sequential differences 
and the method of continued fraction: see interpretation in [2], Sec. 8.2). 
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We remind that the sequential differences Aj of natural numbers Fj and fi are 



Fi+i = max[fi, Ai] , = min[fi, A^] , 



The determinant of the binary quadratic form 



(2) 



_ c^Xf + 2aiXiYi = 1 



is invariant: 



(3) 



biCi + = A. 



For this method (which is algorithmic) one can obtain the expression rj = 6j — Cj + 2aj 
from (2) at the step i. For > 1 the substitutions are: Xj = l^ + Xj+i, = l^+i; for 
Tj < the substitutions are: Yi — X^+Yi^i , Xi = Xj+i (the "ultimate form" becomes 
an identity for X^ — 1, Y^ — 1). Then, substitutions in the reverse sequence give a 
solution of (1). 

We remind the method of continued fraction [8]. Let A > be any radicand 
with \/A =< do;di,d2,---,dL > and L = L(\/A), where a/A =< do;di,...,dL > is 
the continued fraction expansion of VA with period length L. If L is even, then all 
positive solutions of (1) are given by x — Bls-i, s > 1. If L is odd, then all positive 
solutions of (1) are given by x = B2LS-1, s > 1. The sequence {-Bn} is the well known 
recursive sequence given by -B_2 = 1, -B-i = 0,Bn = dnBn-i + -Bn-2 for n > 0. 

The following identity will be used also: 



We remind the well-known statements [2] which prove the procedure of division (for 
the example of + 6^): 

(a) the product of two numbers, each of which is a sum of two squares, is a sum 
of two squares as well (see (4)); 

(b) if the number, which is a sum of two squares, is divisible by a prime which is 
a sum of two squares, then the quotient is also a sum of two squares; 

(c) if the number, which can be written as a sum of two squares, is divisible by 
a number which is not a sum of two squares, then the quotient has a factor which is 
not a sum of two squares; 

(d) if a and b are relatively prime, then every factor of + 6^ is a sum of two 
squares. 

Pell's equation may seem to be fully studied and no new phenomena can be found. 
However, observing the minimal nontrivial solutions of Eq.(l), we see very great 
differences in values of Xi with A^. The answer to the following question is unknown: 
what cases of Ai are the most difficult ones (with maximal values of ajj)? A maximum 
of solutions of Eq.(l) for A between any two nearest squares will be named the local 
maximum, i.e. for k'^ < Al < {k + 1)^ and i — 1,2k we have: ^^[^^^.(A;) > x{Al.). The 
absolute maxima are the such solutions of Eq.(l), that we have x^^^{A'^) > x{Ai) 



(4) 



{y^ - Ax''){Y'' - AX^) = {yY ± AxXf - A{yX ± xYf. 
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for all Ai < A'^. Observations demonstrate, that all local maxima are the cases of 
some prime A = pi and of some double prime A = 2pj, and all absolute maxima 
(for A > 46) are the cases of prime A = An -\- 1 (the record-holders are tested from 
Internet for A < 971853031560). Thus, the study of particular cases of prime A is 
interesting. Since factorization of the double prime corresponding to local maxima 
does not present in the full cycle, these double prime can be named " quasiprime" (we 
do not have the form (2) with = 0, but have one of the characteristic forms for 
prime numbers of A) . Some cases of complex A correspond to one of the characteristic 
forms of prime A (i.e. the case of prime A is rather general one). It is interesting to 
find methods for lowering of volume of calculations (in comparison with well knowm 
methods) for Pell's equation (1) or for a similar equation with arbitrary right side 
(t^ 1), and to obtain explicit formulas. 

We shall use a parametric representation of minimal natural solutions. The pa- 
rameters / and m will denote the basic ones. Knowing the values of I and m (positive 
or zero), we can obtain the minimal positive solution xq for a given number Aq (see 
Propositions 1,4,7, for example). Note that the values of parameters I and m are 
considerably smaller than the value of the solution xq, i.e. we operate with compara- 
tively small numbers. To find these parameters, it is suffice to use a smaller number 
of (initial only) steps than that in (full cycle) the method of sequential differences, 
or, than that in the continued fraction method (for example, see Section 2). The 
parameter i will denote an additional one {i — ±1,±2,...). This parameter deter- 
mines some new natural solution Xi (linearly proportional to i) for some new number 
Ai. Relations (linear in i and Xq), which are determined by the parameter i for fixed 
parameters / and m, will be named "horizontal" relations (see Propositions 2, 5, 8, 
for example). "Vertical" relations are those ones between solutions, for which at least 
one of basic parameters {I or m) is not fixed (these relations are nonlinear in / and 
m). The method of hidden parameters is introduced to find "vertical" relations (in 
the general case for each special case). The method (to find the parameters / and m 
and solutions Xi for Ai) consists in the following: first, we insert (instead of initial pa- 
rameters) an additional (greater) number of parameters; then, resolving some system 
of equations, we obtain a solution (in terms of new parameters); finally, reducing the 
problem to the initial number of parameters, the minimal natural solutions Xi can be 
obtained (as an identity) for Ai. The method is demonstrated in detail to derive the 
Proposition 3 (and is used for Propositions 6, 9). The term "hidden parameters" is 
suggested, since obtaining of the solution (in this procedure) seemed to be impossible 
at a glance (since we seemingly rename parameters only). It is obvious, that the 
relation between a solution (xq, Aq) and other solutions {xi, Ai) does not connected 
with composition of forms in cycles (since the determinants are different for these 
cases). The general solutions demonstrate the relationship between Pell's equation 
(1) and the indeterminate linear equation of type ax — by = 1. 

The purposes of this article are as follows: 

- the method for considerable decreasing the number of calculation steps is demon- 
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strated for Pell's equation (for example, each complicated Fermat case can be solved 
with the help of "a pen and a paper" within a few minutes); 

- by obtaining the solution xq for one can find the infinite number of minimal 
natural solutions Xi for Ai (for filling out the table of the minimal natural solutions 
[Ai, Xi)); 

- the primary aim of this work is to establish some identities for Pell's equation in 
an explicit form. Substituting specific values of parameters l,m,i into these identities, 
the solution {Ai,Xi) of (1) can always be obtained. 

In Section 2 the case of prime numbers A = 4:N+ 1 is investigated. It presents the 
method of decreasing calculations. The infinite number of "horizontal" relations in 
Pell's equation are obtained for Ai — af + bf. The minimahty of solutions is proved. 
The basic ideas, how to obtain the formulae of "vertical" relations in the case of 
A = a'^ + b'^, and the method of hidden parameters are considered in detail. In Section 
3 the case of prime numbers A = 8N + 3 (and the method of decreasing calculations) 
is considered . The connections between {Ao,xo) and {Ai,Xi) are found. Parametric 
solutions are derived. The case of prime numbers A — SN + 7 is studied in Section 
4. In Section 5 the other cases of Pell's equation are discussed. Some connections 
between solutions are demonstrated for each of these cases. The algorithm of "inverse 
calculations" is proposed for the sequential differences. In Section 6 the equation 

- Ax'^ — —3 is considered. Section 7 contains conclusions. 

2 THE CASE OF PRIME A = 4iV + 1 

The case of such prime is of great interest. Notice, that for all special cases, which 
were proposed by Fermat (namely A = 61; 149; 109; 433), the prime numbers are 
4N + 1. Besides, conclusions from this Section can be applied to the case of some 
complex A — a'^ + b'^. 

Proposition 1. Let A be a prime number A = AN + 1 = + 6^ (6 denotes 
an odd number), then the minimal natural solution of Pell's equation (1) is 

(5) x = 2|2Wm + a(Z2-m2)|(Z2 + m2), 
where I and m are taken from the "distinctive form" 

(6) b{f - m^) - 2alm = ±1. 

Proof. The solution of Eq.(l) for a prime number A — AN + 1 is an even number 
(see congruence modulo 4) 

(7) X = 2^SQ 

with coprime odd numbers Q. Since prime A cannot be factorized, two possibilities 
follow: 

{i) y-\^2^^-iAS\ y + l = 2^g% 
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where j = 1 or j = 2T — 1. Subtracting the first equation from the second one, we 
obtain that 2^^^Q'^ ~ 2'^^~^~^AS'^ + 1. This result contradicts the minimal nontrivial 
natural solution (7). 

{ii) y-l^ 2^Q2 ^ ^ + 1 = 2''^-^ AS'' . 

Subtracting the first equation from the second one, it follows, that 2-^~^(5^ = 
22r-j-i^^2 _ rpj^-g equation is impossible to solve modulo 4 for j = 1. Therefore, 
as the result (j = 2T — 1), we obtain 

(8) (2^-ig)2 + 1 = AS'. 

Since the number in the left-hand side of Eq.(8) has the a' + h' type, the number in 
the right-hand side must be of identical type (there exists the procedure of division). 
The value A can be represented as a? + h'. Since one of squares in the left-hand side 
of (8) is unity, the number 5* is a number from primitive Pythagorean triples, that 
is, S — I' -\- m', where I and m are coprime numbers (one of which is odd and one 
is even). The primitive Pythagorean triple is S' — {P — m')' -\- {2lm)'. Substituting 
this into (8), we have 

(a' + b'M' - m'f + {2lmf] = {2^-'Q)' + 1. 

Using the identity (4), we obtain 

b{l' - m') - 2alm = ±1 , 2hlm + a{f - m') = 2^-^Q. 

The net result follows from expression (7). 

Using the expression y = 2A{1' + m')' — 1 and the inequality / > m, we sec, 

that / < \J {y + 1)/{2A) < {\/~Ax + 2)/(2yl). Therefore, the parameters / and m are 
comparatively small ones. We shall demonstrate considerable lowering of calculations 
for the case of ^4 = 61 = 6^ -|- 5^. For comparison, we consider two methods: (I) the 
method of sequential differences; (II) the continued fraction. 

(I) The distinctive form (6) is 5(/^ — m') — 12/m = ±1. The substitutions are: 
1) l = m + C ^ 5C'-12m'-2mC = ±1; 2) C = m + R^ 5R'-9m' + 8mR = ±1; 
3) m = R+D 4:R'-9D'-10RD = ±1; 4) i? = D + V => AV'-lbD'-2DV = ±1; 
5) V = D+E =^ 4:E'-13D'+5DE = ±1; 6) E ^ D+Y^ =^ AY^-SD'+UDY^ = ±1; 
7) D = Yn + G ^ 15Y^-3G' + 8YnG = ±1; 8) G = K + Y^ 20Y^ -3K' + 2KYn = 
±1;9) K = H + Yn^ 19F„2 _ SH' - AHY^ = ±1; 10) if = F„ + X„ ^ 12F„2 _ sx^ - 
IQXnYn = ±1. Therefore, we have = Y^ = 1, H = 2, K = 3, ...,m = 21, ^ 
58. Note that the full cycle for A — 61 contains 64 steps (and the appropriate number 
of reverse substitutions). Substituting parameters m and I into Eq.(5), one obtains 
the well known solution x — 226153980 (instead of usual 64 steps we made only 
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10). In a similar manner, the case of A — 149 can be solved in 12 steps, 18 steps are 
necessary for A = 109, and so on. The question of representation of prime A = 4N+1 
as + 6^ can be bypassed (sec Section 5 below: for "inverse calculations" we can 
start from the form 3y^ — 12x^ — lOxy = 1). 

(II) Now we consider the method of continued fraction. The continued fraction 
expansion of y/A gives the following expression at the step i: 

(9) l<i<L + l. 

Note that these coefficients determine some form (2). In this case, we have: 

= { — ly'^^al, 62?+! — — T2i-i C2i+2 — Cli+1 — T^2i+1- 

The method for lowering of calculations consists in the following: we seek two ex- 
pressions in the continued fraction expansion with rj = rj_i. Therefore, we have 
a — Oj, b — Tj. The parameters / and m is / = -Bj-i, m = Bj^2, where {B^} is the 
recursive sequence given by B_2 — 1, -B-i — 0,Bn — dnB^-i + -Bn-2- Thus, 

rj = Tj-i =^ a = a*, h — rj, I — Bj_i, m = Bj_2. 

We demonstrate the case of A = 61. 

v^=7+(V61-7) = 7+ ^ , .0 = 1 + -^ — = 1 + 



61 + 7 12 12 V61 + 5 

V6T + 5 ^ I - 7 ^ 4 761+7 ^ =3 I 9 



3 3 ^61 + 7' 4 4 v^ + 5' 



61 + 5 V61-4 5 \/61 + 4 V61 - 6 5 
= 1 + ^ = 1 + ^=^ — -, = 2 + = 2 + 



9 9 x/6T + 4' 5 5 ^ + 6' 

We see that 6 = 5, a = 6, = 1, Si = 1, S2 = 5, ^3 = 16, m = = 21, 1 = B5 ^ 58. 
Substituting 1,171, a, b into Eq.(5), the minimal positive solution x follows. Instead 
of 12 steps for the continued fraction expansion of v^6T we use 6 steps only, and 
instead of calculation of B21 (it is the solution - see Introduction), we calculate B^ 
only. Therefore, the number of calculation steps can be considerably decreased. 

Proposition 2. Let xq be the minimal natural solution of Eq.(l) for Aq — al + b^ 

and definite parameters l,m from Eq.(6); then xi are the minimal natural solutions 
of Pell's equation for Ai andi ^ —1, where 

(10) Xi^\ xo + 2iS^ \, S^P + rri^, 

(11) A^Ao + i^ + i'S'. 
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Proof. In fact, parameters I and m are fixed ("horizontal" relations). Introducing 
substitutions 

(12) k = f-m^>0, t = 2lm, 
we rewrite (6) as 

(6') bk-at^±l. 
It follows from this: 

(13) bi — bo + it, Qi — ao + ik, i = ±l,±2, ... 

Substituting (13) into expression (5), we obtain Xj = 2|6ji + aj/c|(/^ + m^). Taking into 
account substitution (12) and the identity k"^ + 1"^ — S'^, one gets the net formulae 
(10) and (11). The values {l,m) are determined by {k,t). It is well known that there 
exists one minimal positive solution (A;o,to) of Eq. (6') with |A;o| < |a|, |to| < \b\ for 
the given a and b. Therefore, the same values {ko, to) are minimal ones for and bi, 
i ^ —1 (see expression (13)). The solution minimality is proved. 

Note, that if Xq is some nonprimitive solution, then x_i can be a primitive so- 
lution (if |a_i| > |ao| and > \bo\). For example, we can find from the case 
of Aq = 2 (parameters / = 2,m = 1 give the non-minimal solution xq = 70): 
Xi=\70 + 250i I for = 2 + Ui + 25f. The solutions arc (A„ xf, i): (13,180:-!) , 
(41, 320; 1), (74, 430; -2), (130, 570; 2), etc. Now we describe the general method for 
finding "vertical" relations. 

The method of hidden parameters. Wc introduce additional parameters (that 
is, the problem must be complicated at first), and then the relations between these 
new parameters must be applied. For the case oi A = + b"^ we suppose, that all 
four values 1,171, a and b are unknown, and shall seek the solution of condition (6) as 

(14) l = gn + r, m = dn + h, b = pn + z, a = qn + u, 

with new unknown integer parameters n,g,r,d,h,p,z,q and u. Seemingly, the sub- 
stitution n = into final identities is to no avail, because it implies only rewriting 
of symbols, which it is not the case. This is the reason, why the notion of "hidden 
parameters" appears in the title. 

We suppose that n can be arbitrary in an ultimate solution and group Eq.(6) in 
terms of n: 

(15) (n^) : g^p-d^p-2gdq = 0, 

(16) (n^) : 2grp - 2dhp + g'^z - d^z - 2ghq - 2drq - 2ugd = 0, 

(17) (n^) : pr^ - ph'^ + 2grz - 2dhz - 2rhq - 2ugh - 2dru = 0, 



7 



(18) : z{r^ - h'^) - 2urh = ± 1. 

Our aim is to find any solution of Eq.(6), namely, a and b for distinct parameters I 
and m. Then, using "horizontal" relations (13), we shall obtain all solutions {Ai,Xi) 
of Pell's equation. One can put in (15): 

(19) P = gdv, q = ^v{g^ -(f), 

where v is an arbitrary unknown integer. Substituting p and q into (16), (17) and 
resolving the linear system of equations in terms of v and u, one gets: 

g^ + (P z\h(f - rg^ - 2>gdird - hg)] 

^ ^ hg^ + zd^ 2{hg^ + rd^) 

It follows from Eq.(18), that z and u are relatively prime integers; z is an odd value 

(21) z^hg^ + rd^. 
Substituting z in Eq.(18), one obtains the condition 

(22) rd-hg^ ±1. 
Thus, the final expressions are as follows: 

77 1 

(23) 6o = ngd{g^ + d^) + hg' + ro^^ ao = -(^^ - d^) - -[hd' - rg' T 3^ci] 

with a single condition (22). If for arbitrary coprime numbers r and h we shall know 
any g and d from condition (22), then the numbers oq and 6o can be calculated from 
expression (23); therefore, taking into account the "horizontal" relations (13), all 
solutions of Pell's equation with parameters / and m are found. 

The remark on "hidden parameters" is in order. For n = we have the following 
solution (which can be proved by direct substitution) from expression (23) despite 
the fact, that it seemed to be impossible at a glance (we rename r ^ Z, /i — > m). 

Proposition 3. The minimal natural solutions Xi of (1) for different — of + hi 
can be taken from expressions (5), (10), (11) with 

(24) bo = mg' + ld', gq = ~[md' - Ig' =F 3gd], 

(25) Id-mg^ ±1. 



Selecting the identities, which represent the unity in various forms (and comparing 
with (25)), and using formulae ((5), (10), (11), (24)) derived above, the representation 
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of Pell's equation in form of identities can be written, and for specific parameters these 
identities present the solutions which arc of interest to us: 

(a) Putting d = 1, we have (see (25)) / = mg ± 1 (for the odd number m the 
number g is odd; for the even m the value g is arbitrary), 



[2(m^ + m'^g'^ + 1 ± 2mg)2m{mg ± l)ai + {m'^g'^ + 1 ± 2mg - m^)6i]^ + 1 = 



For example, we present (from here) some minimal natural solutions for A < 150 (see 



for g = 3, m = 2, i = —2 and at the top signs we have A = 89, x = 53000; 

if g = 5, m = 2, i = —6 (the top signs), then x = 267000 for A = 73; 

for g( = 3, m = 3, i = — 2 (the bottom signs) we obtain A = 113, x = 113296. 

For other identities, which can also be written with the help of "a pen and a 
paper" only (with using formulae (5), (10), (11), (24), (25)), we present parameters I 
and m, A and x: 

(b) Putting in (25) d = g — 1 (the bottom sign), it turns out that some identity 
can be written with / = g{2T + l) + 1, m = g{2T+ 1) — 2T; for T = 3, gf = 3 we have 
Z = 10, m — 7, and with "a shift" i — —2 the minimal natural solution for A — 137 
is X = 519712. 

(c) Let g = rd^l , m = 2dT + 1,1 — rm ^ 2T (to obtain some identity in (25)); 
for T = 3, d = 2, r = 2 (top sign) the parameters are m = 13, I — 20, and with "the 
shift" i = -lwe find A = 97, x = 6377352. 

(d) We substitute g = J + T( J- 1), d=J-l, r = TJ+J + 1, l = 2gni+r, m = 
2dni + J (we have an identity from (25)), and for J = —3, T = 2, m = — 3 the 
parameters are m = 21, / = 58; the latters (with "the shift" i = 13) describe the case 
A = 61, X = 226153980. 

(e) We have / = rii + Tm, g = 1 + dT, m = riid^l ((Eq.(25) becomes an identity); 
it follows for Til = 5, d = 6, T = 3 (the top sign), that m = 29, I = 92; the latter 
parameters with i = —41 satisfy the case of A = 149, x = 2113761020, etc. 



bi = mg^ + mg ± 1 + i{2m^g ± 2m), 

ai = —-\rn — mg ^g ^3g\ + i{mg +l±2mg — m), 
and the three-parametric identity (for arbitrary values of m, g, i) is 

{mg^+mg±l+2m'^gi±2mif+ ^) ^ 9{9 j^m'^g'^i±2mgi+i-imA x 




[2{af + hDirn'g'' + 1 ± 2m^ + rn'f - if. 



formulae (5), (10), (11), (24), (25)): 
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3 THE CASE OF PRIME A = SN + 3 



Proposition 4. Let A be a prime number A = 8A^ + 3 = + then the minimal 
natural solution of Pell 's equation (1) is 

(26) X = I Ahlm + a\f - 2m^\ \ ■ {f + 2m^), 
with I and m taken from the distinctive form 

(27) h\f - 2m^\ - 2alm = ±1. 

Proof. The minimal nontrivial solution is an odd number in this case (see the case 
of A = 4Ar + 1) 

(28) X = QS. 

Using the simplicity of A, two possibilities follow again: 

(i) y-l^AS\ y + l^Q^ 
Subtracting the first equation from the second one, it turns out that 

(29) Q^^AS^ + 2. 
However, this equation is contradictory to solve modulo 8. 

(ii) y + l^AS\ y-l^Q'. 

(30) = + 2 • 

Since the number in the right-hand side is representable as a sum of a square and 
a double square, the number in the left-hand side has similar representation (the 
procedure of division). We present A as + 26^. Since one of squares is unity, S may 
be represented as S = P + 2m^ and S"^ — — 2m^Y -\- 2(2Zm)^. Using expressions 
(4), (30), the proposition is proved. 



Note, that this case can be apphed to some cases of complex A — -\- 2h'^. 
Using the expression y — A{1'^ -\- 2m?Y — 1, we see, that l,m < {y + 1)/^ < 

{\fAx + 2) I A. We have considerable decreasing of calculations again (for the method 
of sequential differences), since we use some part of the cycle only and calculate 
comparatively small parameters Z, m. The question of representation of ^4 as + 2\? 
and of choosing the sign (for removing of the modulus) can be bypassed with using 
the inverse calculations (see Section 5 below). 

We present the method for decreasing of calculations for the method of contin- 
ued fractions. Using Proposition 4 and the correspondence of the form (2) and the 
expression (9), we obtain the following result. 
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(I) If we have 



Vj = 2rj_i ^ a = Gj, b = rj^i, m = Bj_2, I = Bj_i. 
(II) If we have 

Tj-i — 2rj =^ a — Up b — rj, m — Bj^i, I — Bj_2- 

Everywhere {-Bn} is the recursive sequence given by i?_2 = 1, -B-i = 0, = dnBn-i + 
Bn-2- We demonstrate the case oi A~ 139. 



^l39+ll' 18 ^139 + 7' 5 ^139 + 8' 

v^ + 8 ^ 6 + 7 ^ 3 
1 + ^^ , - — = 3 + 



15 ^139 + 7' 6 ^139 + 11' 

Wc sec that (II) 6 = 3,a = ll;Bo = = 1,-02 = 4,/ = S3 = 5,m = E4 = 19. 
Substituting 1,171, a, b into (26), the minimal positive solution x = 6578829 follows. 
Instead of 19 steps for the continued fraction expansion of yl39 we use 5 steps only, 
and instead of calculation of Bn (it is the solution - see Introduction), we calculate 
B4 only. Therefore, the number of calculation steps can be considerably decreased. 



Proposition 5. Let xq be the minimal natural solution of Eq. (1) for Aq = a\ + 26q 
and definite parameters l,m from condition (27), then Xi are the minimal natural 
solutions of Eq.(l) for Ai andi ^ —I, where 

(31) Xi = \xo + iS^\, 

(32) Ai^Ao + 2i^+i^S\ 

o 

Proof. Designating k = — 2m?\, t = 2lm, the condition (27) can be rewritten 
as (6'). Because Eq.(6') is satisfied by Gq and bg, all solutions are (13). The solution 
can be written as Xi — \aik + 2bit\S, A^ — af + 2bf. Because k'^ + 2t^ — S'^, the net 
formulae follow. The proof of solution minimality is vahd (see Section 2). 



For example, (using the nonminimal solution Xq = 12 and parameters 
/ = l,m = 1) we can obtain the minimal natural solutions Xi = |12 + 
27i| for Ai = 2 + 8i + 9i^. The {Ai,Xi;t) solutions (from (31), (32)) are: 
(19, 39; 1), (54, 66; 2), (107, 93; 3), (22, 42; -2), (59, 69; -3), (114, 96; -4), etc. 

Proposition 6. The solution of Eq.(l) for the case 0/ A = 8Ar + 3 = aj + 2bl is 
presented by expressions (26), (31), (32) with 

(33) 60 = ^mgl + Id^ , ao = \md^ - 2lgf ^ 3gid\ ; 
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(34) ld-2mgi^ ± 1. 

Proof. The direct substitutions prove this proposition. However, we shall outline 

the derivation of the statement with using the method of hidden parameters. Sub- 
stituting (14) into condition (27) (we substitute sign = ± instead of |...|) and 
grouping power series in n, we have four equations again. Substituting 

(35) p = sign-gdv, q = ^v{g'^ - 2d'^), 

and resolving the linear set of simultaneous equations in terms of v and u, we find: 

. g'^ + 2d'^ . z[Ahd^ - rg^ + Qgdihg - rd)] 

(36) V = sign ■ z , - sign 



hg^ + 2rd^ ' ^ 2{hg^ + 2rd^) 

In order that Eq. (27) be solvable, it turns out that 

(37) ^= l(/,^3^2rrf3), 
where g — 2gi. The final condition is 

(38) rd — 2hgi = =F " 1 
The solution is: 

(39) ao = sign ■ [{Ag'l — d'^)n — [hd^ — 2rgl ± sign ■ 3gid]], 

(40) bo = 2gid{2gl + d^)n + Ahgl + rd^. 

"Horizontal" relations (31), (32) provide all solutions for the given parameters. Al- 
though we have sought again only some of solutions, all solutions of Eq.(l) can be 
obtained with the help of condition (38) in this case. Letting again n = and rename 
r — > /, /i — > m, we have the net result from the proposition. 



For example, some identities for (1) can be obtained in the following way. 

(a) II d — 1, then (from (34)) I — 2gim ± 1; for gi — 2, m — 1 and the top sign 
one gets / = 5; for "the shift" i — —A we have (from (26), (31), (32), (33)) the case of 
A = 67, x = 5967. 

(b) If gii = 3, m = 1 (the top sign), then / = 7, and with i = —8 we have the case 
of ^ = 118, = 28254. 

(c) We substitute in (34): d = 2Tgi ± 1, and Z = 2giK + 1, m = Tl ± K. Let 
gi = 2, K = 1 m the identity; then / = 5, and for T = 4 (the bottom sign) the 
parameter m is: m = 19; with "the shift" i = —92 the solution follows: x = 6578829 
for A = 139. 
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4 THE CASE OF PRIME A = SN + 7 



On one hand, the solution x cannot be an even number (see Section 2). On the 
other hand, Eq.(30) is intractable to solve modulo S ior A — SN + 7. Therefore, we 
obtained the solution (28) with coprime odd parameters which can be found from the 
condition (29) by sequential differences. We note, that the representation of prime 
A = 8N + 7 as a^ — 26^ follows from (29) (the existence of the solution and the division 
procedure). Substituting A = — 2b^ and S = P — 2w? into (29); thus the following 
proposition is proved: 

Proposition 7. The minimal natural solution of Eq. (1) for prime ( and some 
complex) A — — 2b^ is 

(41) X = I (0(^2 + 2m^) - mm) {f - 2m^) \ , 

(42) 2alm - b{l^ + 2m'^) = ±1. 



Since / > m (wc substitute m = l — z, further I = m + z), the distinctive form (the 
sum of coefficients at the squares of parameters equals the coefficient at the product 
of parameters) is 

(43) {2a-3b)m^ -bz'^ + 2{a-b)mz^±l. 

Writing y = A{1'^ - 2m'^f + 1, wc obtain /, m < ^jy/A < ^ {VAx + 1)/A (i.e. the 
parameters / and m are comparatively small again). Using (2), (9) and Proposition 
7, the method for decreasing of calculations for the method of continued fractions can 
be suggested. If we have 

rj_i + Tj = 2a* =^ b = rj_i, a = a* + r^.i, m = -Bj-2, I = Bj_i + Bj_2 

with the recursive sequence {-B„} given by B_2 = Ij B_i = 0, i?„ = dnBn-i + -Bn-2- 

One can write the relations between {xo,Af)) and {xi,Ai), so that all substitu- 
tions in the method of sequential differences are identical. Designating t = 2lm, k = 
P + 2m^, the condition (42) can be rewritten as (6'). Because Eq.(6') is satisfied 
by ao and 60 ; all solutions can be taken from expression (13). The solution can be 
written as Xi — \ {aik — 2bit)S |, A^ — — 2bf. Because S'^ + 2t^ — A;^, we proved 
the following proposition: 

Proposition 8 All "horizontal" relations between minimal natural solutions xq 
for Aq and Xi of Pell's equation (1) (and i ^ —1 ) for the case of Ai — — 26| are 
presented by the formulae 

(44) Xi = \xQ + iS'^\, S = f- 2m^ 
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(45) Ai^Ao + 2i^+fS\ 



For example (we use parameters / = 3, m = 1 for nonminimal solution x = 70, A ~ 
2),Xi^ |70 + 343i|,A = 2 + 20i + 49i2, i.e. {Ai,Xi,i) : (31, 273; -1), (71, 413; 1), etc. 

All "vertical" relations can be found by the method of hidden parameters. We 
present the results for expression (14). 

Ahd^ - rg^ ± 6gd 2d^r - hg^ 



u 
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For n = 0, g = 2gi (with r ^ l,h ^ m) we have the following proposition, which 
can be checked by the direct substitutions. 

Proposition 9. The solution of Pell's equation (1) for Ai — — 26? is given 
by expressions (41),(44))(45) with 

(46) ao = md^ - 2lgf ± 3gid, bo = Id^ - ^mgf, 

(47) Id - 2mgi = ±1. 



We present the examples of identities below. 

(a) The condition (47) is an identity for d = 1, m = 1, Z = 2g'i + l. For gi — 3,i — 7 

the case A = 103, x = 22419 follows from (41),(44)-(46). 

(b) The condition (47) becomes an identity for = l,m — 2gi,l — Ag\ ± 1. If 
g^ = 2,i = l, then (see (41),(44)-(46)) A = 127, x = 419775. 

5 OTHER CASES 

The rest case of ^4 = P1P2, where pi and p2 are coprime numbers (however, they can 
be complex numbers), can be broken down into the following ones. 

(A) The solution a; of (1) is an even number x — 2QS. It follows from 

y+1^2p,S^ , y-l = 2p2g', 

that the distinctive form is 

(48) Pis'" - p2Q^ ^ 1. 

For the method of continued fractions, this form does not present in (9). However, 
expression some parameter in terms of other one leads to the form, which is presented 
in (9). For the decreasing of calculations, if we have 

Krj^a], X = l,2,..., =^ Q ^ Bj^^, S ^ KBj^^ + Bj_2: 
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where the recursive sequence is B_2 = 1) = 0, i?„ = dnB^-i + -Bn-2- 

Since p\ — poi + iQ^, P2 — P02 + iS^, these "horizontal" relations give 

2 

(49) Xi = 2g5, A^Ao + i^^ + i{poiS^+po2Q^). 

The following "vertical" relations can be checked by direct substitutions: 

(50) poi — m^{3Ql — Sm), po2 — l^i'^Sm — Ql), 
with the condition 

(51) Ql-Sm^ 1. 

The solution always exists for coprime numbers S and Q; that is, we have obtained 
all relations for Xi and = P01P02 (see (49), (50)). 

(B) The solution of (1) is an odd number x — SQ. The condition is 

PiS^ - = ±2. 

Let pi > p2, Q = S + 2q, and one obtains the distinctive form with equal coefficients 
at one of squares and at the product of parameters 

(52) PLZPls' - 2p2q' - 2p2Sq = ±1. 

For the decreasing of calculations in the method of continued fractions, we seek the 
following condition in (9): 

rj = 2a* X ^ Bj_i{Bj^i + 2Bj_2) 

with the well known recursive sequence {Bn}- 
The "horizontal" relations 

(53) Xi = SQ, Ai^Ao + i{pmS'^ + Po2Q^) + i^xl 
with the "vertical" relations 

(54) poi ^2m^{2,Ql- mS) , po2 = {2,mS - Ql), 

(55) IQ-mS = 1 

present all solutions of Pell's equation for this case (and can be proved by direct 
substitutions) . 

We write the binary quadratic form (2) with 6„ > 0, c„ > and 

(56) bn-Cn + 2an = 1 
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(the sequential differences lead to — — 1 for some step n). The net formulae 

are 

(57) A-k{k+l). 

(58) bn^{k + If - A, 

(59) c„^A- k'^. 
It follows from 6„ > 0, c„ > 0, that 

(60) k^ <A<{k + If. 

This condition for the given A (not a square) uniquely determines the value k and, 
therefore, the values a„, c„. It can be proved that (2) can be reduced to (1). 

Inverse calculations (for the method of sequential differences). We believe that 
the "coordinates" of a given number A are known (the nearest square k"^ and the 
"distance" from it). Therefore, using (57)-(60), one can write the "ultimate form", 
which represents the "one" (by substitutions X„ = 1,1^ = 1). Using the symmetry 
of cycle with respect to the forms with +aj and — Oj, we change the sign of a„ and 
obtain the "start form". By the sequential differences method we obtain the "final 
form" . This is one of five types of the form discussed above: 

(I) the coefficients at the squares of parameters are equal (see (6)); 

(II) the coefficients at the squares differ twice from each other (see (27)); 

(III) the sum of coefficients at the squares of parameters equals the coefficient at 
the product of parameters (see (43)); 

(IV) the product of parameters is missing from the form (see (48)); 

(V) the coefficient at one of squares equals the coefficient at the product of pa- 
rameters (see (52)). 

Changing the sign of the product of parameters, we obtain the "distinctive form" 
(with real sign), that is, the real sign in the solution follows. To find the "distinctive 
parameters" , it is necessary: 

1. to find the "start parameters" Xn and Yn in terms of the "final parameters"; 

2. to express the "final parameters" in terms of the "start parameters"; 

3. to change the sign at one of the "start parameters" (either X„ or Y.^) in both 
expressions, to substitute X„ = F„ = 1 and to take the magnitudes. As a result, we 
obtain the "distinctive parameters" (and the solution of (1) for the given number A). 
Notice, that beginning from the "start form" we can also found the representation of 
number A as + 6^, or as + 26^, or the expansion of A into two coprime factors 
(the same result can be obtained in continued fraction expansion). 

For example, we outline the case oi A = 103. The "start form" (10^ < 103 < 
112, A; = 10) ^gy^2 _ 2,xl + 14X„y; = 1. We have the "final form" (III) at the 5th 
step: 13mf — ?>zl — IdmiZi — 1. We obtain 6 = 3, a — 11, Xn — mi, Yn — 5Xn + z\. 
Therefore, z = 6, m = 1, Z = 7 and it follows from (41): x = 22419. 



16 



6 SOME RELATED PROBLEMS 

The methods described above permit to exphcitly solve some related problems in the 
general case (for example, y'^ — Ax'^ — ±j^; ±2j^, etc.). We consider the case (see [7]): 

(61) - Ax^ = -3. 

Since there exists the procedure of division for the numbers + 3j^, two possibihties 
follow for odd A (the case of AA can also be reduced to them). 

(i) The minimal natural solution a; of (61) is an odd number for A — + 3b^ — 
AN + 3 (with an odd number h) in the first case: 

(62) x = ^^ + 3m^ y = I (a|Z^ - 3m^| - 66Zm) I, 
with the condition 

(63) 2alm-h\f -?>m^\^±l. 

To decrease calculations for the method of continued fractions, we obtain the following 
result. 

(I) If we have 

Tj — 3rj_i ^ I — Bj_i, m — Bj_2, a = a!ij, b — rj_i; 

(II) if we have 

Tj-i = Srj I — Bj^2, TTi — -Bj-i, b — rj, a — a* 

with the recursive sequence {-B„} : B_2 = 1, B_i = 0, = d^Bn-i + -Bn-2- 
All "horizontal" relations can be taken from 

(64) = ao + — 3m^|, bi — bo + 2ilm, Ai = a^ + 3b^. 

Using the method of hidden parameters, one can find the parametric solution of this 
problem (which can be proved by direct substitutions): 

3 

(65) ao^ -[d^m + 9dmg^ -Sd^lg -3lg% bo = 9g^m + ld^, 

where the condition for d and g is 

(66) dl - 3gm = ±1. 

(ii) The second case with an even solution of (61): x = 2(/^ + 3m^), takes place for 
A = Al + 3Bq = AN + 1 {Bo is an even number). In this case {x 2xi, A AAi) 
all previous formulae are valid with substitutions 

(67) a^Ao± 3Bo, b^A^T Bq. 
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We consider the case A — 1729 = 4A'" + 1 (see [7], for comparison). 

r—r— 27 v/69l6 + 83 ^ 25 
V4 • 1729 = 83 + , , — = 6+ — 



v/69T6 + 83' 27 + 

^6916 + 79 , 75 
= 6 + 



25 ^6916 + 71' 

We see that (I) m = Si = 6, Z = S2 = 37, a; = 2954, y = 122831. We outhne the 
parametric identity which includes this case (see (62), (64), (65)): the condition (66) 
is an identity for d = 1,1 = 3gm + 1 (for m = 6, g = 2,i — —1 one obtains from 
(64), (65): 6 = 25, a = 71, and the case follows from (67)). 



7 CONCLUSIONS 

The general case of Pell's equation (1) can be broken down into the special cases, and 
in each of these cases the number of calculation steps can be considerably lowered for 
the method of sequential differences and for the method of continued fractions (the 
minimal natural solution can be found with the help of "a pen and a paper" within a 
few minutes). The item-by- item examination of particular cases may seem increasing 
the volume of calculations. However, it is not the case. The author computed Pell's 
equation (in Python): (a) with using the method of continued fraction (standard 
algorithm), and (b) with using (in addition) the particular cases and formulas from 
this article. We have the following results. For all A < 1000000 it is necessary 
about 3 hours and 2 minutes for the standard algorithm and 1 h. 43 m. for the 
modified algorithm; for all A < 2000000 it is required about 8 h. 50 m. for the 
standard algorithm and 4 h. 50 m. for the modified algorithm. Thus, the lowering of 
computation time is not less than 1.828 times in the last case (and 1.827 times, if we 
seek and write the absolute maxima also). Since the maxima significantly increase 
with increasing of A, this result can be improved (the basic time surplus we have 
from the most difficult cases). Computer programs (from Internet) use the method of 
continued fraction. The most fast algorithm (in Mathematica) uses additionally one 
particular case only: yf — Ax^ = —1 (factually, Eq.(8) from Section 2, but not Eq.(6)). 
Therefore, the suggested algorithm is more effective in this particular case also. The 
methods, which permit to decrease the scope of calculations for finding the table of 
minimal natural solutions of Pell's equation (1), are investigated in the work. Using 
these methods, it is possible to estabhsh the relations (to deduce the appropriate 
formulae) between a solution xq for Aq and the infinite number of minimal natural 
solutions Xi for Ai in various cases. 

The establishment of "horizontal" (both basic parameters are fixed) and "vertical" 
(basic parameters are variable) relations permits to write a parametric representation 
of Pell's equation (1) as identities. With the help of "the method of hidden param- 
eters" all relations can be found. This method consists in the following procedure: 
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first, it is necessary to complicate the problem by inserting an additional number 
of new parameters; then, by obtaining relations between these new parameters, the 
appropriate solutions can be found; and, finally, by reducing the problem to the ini- 
tial number of parameters, the solutions can be represented as identities by means of 
initial parameters. Here these relations permit to find the minimal natural solution 
of Pell's equation (1) (for the given Aq) and to present Pell's equation (1) as a para- 
metric identity in the explicit form (for infinite natural numbers of {Ai,Xi), including 
{Aq,Xo)). These relations do not connected with composition in cycles, but Pell's 
equation (1) is connected with indeterminate linear equation ax — by = 1. 
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